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We observe a resonance in the conductance of silicon tunneling devices with a 5— doped barrier. 
The position of the resonance indicates that it arises from tunneling through the state of the 
boron atoms of the 5— layer. Since the emitter Fermi level in our devices is a field-independent 
reference energy, we are able to directly observe the diamagnetic shift of the level. This is con- 
trary to the situation in magneto-optical spectroscopy, where the shift is absorbed in the measured 
ionization energy. 



The smallest semiconductor device with potential func- 
tionality is a semiconducting nanostructure with a single 
dopant atom. The properties of such a structure are most 
prominent at low temperature, where the electron or hole 
is localized at the parent donor or acceptor, respectively. 
Manipulation of the wave function of the charge carrier 
at the dopant atom with the electric field of a gate is the 
obvious tool to influence the properties of the nanostruc- 
ture. A good first order description of the properties of 
dopant atoms is given by the hydrogen model^. The re- 
sulting Bohr radius of up to about 10 nm sets the size of a 
dopant atom and thus of a single-dopant-atom functional 
device. This radius is much larger than for the hydrogen 
atom, due to scaling for the effective mass of the carrier 
and the dielectic constant of the semiconductor. Such 
dimensions are accessible with e-beam lithography and 
scanning probe techniques. This opens the way for what 
can be called atomic scale electronics inside a semicon- 
ductor. Silicon is very attractive for this purpose because 
of its highly developed fabrication technology. 

A beautiful example of atomic scale electronics inside 
silicon is the quantum computer proposed by Kane^ and 
by Skinner et al^. Both the fabrication and electrical 
operation of the qubits of this computer rely on control 
at the level of individual phosphorus donors. The highly 
controlled dopant engineering required for atomic scale 
electronics inside silicon is being worked on by several 
groups^. Doping at the atomic scale for application 
purposes is thought to be feasible, but development of 
this technique will be time consuming. 

A more direct way to a single dopant atom in silicon, 
albeit less controlled concerning exact positioning, is to 
use (5— doping and conventional nanostructuring. Single- 
dopant-atom structures fabricated in this way will yield 
physics relevant for future devices fabricated with atomic 
scale doping techniques. To some extent this approach 
has been followed already, by including a (5— layer of 
dopant atoms in the well of double barrier diodesS. Here, 
we report a transport study on 5— doped silicon tunneling 
devices grown with a single barrier. The dopant atoms 
induce zero-dimensional atomic quantum wells, giving 



many identical double barrier systems in parallel. This 
device, a precursor of a single-dopant-atom device of a ge- 
ometry close to a qubit of Kane's computer, shows very 
interesting transport properties. In particular, we find 
a conductance resonance due to tunneling through the 
boron impurities of the 5— layer in the barrier. The posi- 
tion of the resonance and its magnetic field dependence 
indicate that it originates from tunneling through the B+ 
state of the boron impurities. So far, this state has only 
been observed in spectroscopic studies using photons or 
phonons and not in an energy resolved transport experi- 
ment like ours. 

We fabricate the (5— doped devices from a layered struc- 
ture of the type p+ Si(500 nm)/p- Si(20 nm)/(5/p- Si(20 
nm)/p+ Si(500 nm). Boron is the dopant, for the layers 
and the (5— spike. The structure is deposited by chem- 
ical vapor deposition in an ASM Epsilon 2000 reactor 
on a Si(OOl) substrate with low doping, using SiH4 and 
B2H6 as precursor gasses. The (5— spike of areal den- 
sity 1.7 X 10^^ cm~^ is centered in the lowly doped 40 
nm thick tunnel barrier. The p"*" layers are degenerately 
doped {Nb — 10^^ cm^'^) and serve as contact layers. 
The devices are square mesas, 100, 200, 300 and 400 
/im at a side. They are dry-etched in an SFg plasma. 
The etch mask is the Al-l%Si top contact of the mesa, 
which is sputter-deposited through a shadow mask. The 
SFg etch is stopped just after the bottom p+ layer has 
been reached. A second shadow mask, aligned with re- 
spect to the mesas, is used for sputter deposition of Al- 
l%Si contacts to the bottom layer. The final step is a 
400 °C anneal in N2/H2 of the Al contacts to the p+ 
Si, using rapid thermal processing. A high device qual- 
ity is apparent from resistance scaling with mesa size. 
At room temperature the resistance is dominated by the 
two-dimensional spreading resistance of the bottom layer 
between the mesa and the Al contact, while at low tem- 
perature it is determined by the barrier in the mesa. 
These devices, which have a metal-insulator-metal struc- 
ture, are the simplest all-silicon tunneling devices. 

We measured the doping profile in the structure with 
secondary ion mass spectroscopy (SIMS). Figure Qfa) 
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Figure 1: (a) SIMS profile of the boron concentration. Zero 
depth is the surface position. For optimum resolution the 
top layer is only 20 nm thick. Contact layers, barrier and 
5— layer are clearly visible, (b) Depicts the profile of the 
valence-band edge, the hole energy increasing in the upward 
direction. Fermi seas of emitter and collector extend up to 
the barrier, of which the thickness is defined with the criterion 
Nb = 4 X 10^* cm~^, the concentration of the metal- insulator 
transition. Symbols are discussed in the text. 



gives the result, in which the 5— layer, the barrier and 
the contact layers are clearly discernible. The (5— layer 
is about 6 nm wide and has a peak concentration of 
Nb = 5 X 10^^ cm"^. The contrast of the (5— layer and 
the background doping {Nb ~ 10"'^^ cm~'^) in the barrier 
is somewhat weak, although for the barrier B2H6 is only 
applied during (5— doping. This is due to boron diffusion 
during the silicon growth, out of the p"*" bottom layer 
and out of the (5— layer. The profile of the valence band 
edge is depicted in Fig. ^b). The tunnel-barrier height 
isi (j)B = AEy — Ep, i.e. the valence-band contribution 
AEy to the band-gap narrowing of emitter and collector 
minus the Fermi energy Ep of these device layers. 

Electrical measurements were performed in a flow cryo- 
stat equipped with a 14 T superconducting magnet and 
in a '^He cryostat. We use standard lock-in techniques 
to measure G — V curves, i.e. curves of the differential 
conductance versus bias. In the G — V curves, at 4.2 K 
and below, a tunneling resonance is present around 10 
mV, superimposed on a dominant background. In Fig. [21 
we show curves of a 400 fim device, for temperatures 
between 0.5 and 12.5 K. The weak resonance at 4.2 K 
becomes a clear peak at lower temperature, which is not 
yet saturated at 0.5 K. At this temperature the FWHM 
of the peak is about 1.5 mV. The overall behavior as 
sketched in Fig. [3 is present for all measured devices, 
which are of different size and come from different fabri- 
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Figure 2: Conductance curves of 400 [xm device 1/888/1, for 
the temperatures listed. With decreasing temperature the 
conductance decreases (for the lower temperatures most clear 
at higher biases) and a resonance at ±10 mV develops. Be- 
tween 1.5 and 0.5 K the resonance still grows. 



cation runs. The peak is absent for devices that have no 
(5— layer but are otherwise identical. This indicates that 
the peak originates from the boron atoms of the 5— layer. 

In the inset of Fig.|31we plot the response of the spectra 
(background subtracted) to a magnetic field, for the same 
device as in Fig. |2 for one bias polarity* and for fields 
up to 14 T oriented perpendicular to the layers. With 
increasing field the peak shifts to higher bias and becomes 
broader and lower. In the main panel of Fig. |31 data 
points of the level shift A£' = ieAVros deduced from the 
resonance shift AVrcs (see below for the relation between 
level and resonance positions) follow a weak quadratic 
function of the field, the shift at 14 T being about 1 
mV. This behavior does not depend on the orientation of 
the field with respect to the crystallographic axes of the 
device. 

The conductance peak is attributed to resonant tun- 
neling through the state of boron impurities in the 
(5— layer, each of which provides for a conductance chan- 
nel. The B"^ state is an acceptor counterpart of the more 
generally known state^. It forms when a second hole 
is weakly bound to a neutral acceptor, in our case neutral 
boron B°. In zero magnetic field the B+ state is a sin- 
glet state (as the D~ state is), which is analogous to the 
negative hydrogen ionS (H~ ion). The separation of the 
B" ground-state level and the B"*" level results from the 
Coulomb interaction between the holes. In our devices 
the B° ground state is deep below the Fermi level, so 
that it is permanently occupied. Therefore, higher levels 
of the B° single hole spectrum are not available as step- 
ping stones for holes tunneling from emitter to collector. 
This means that the B+ state is the only candidate which 
can induce the resonance. This is unlike an optical ab- 
sorption experiment in which higher levels may always 
come into play via transitions from the ground statei. 

/^From photoconductance spectroscopy on Si samples 
with low doping level it is known^" that in zero magnetic 
field the binding energy E^(^) of the extra hole on an 
isolated B+ ion, i.e. the amount of energy required to 
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Figure 3: Magnetic field induced shift of tfie discrete energy 
level in the barrier of device 1/888/1, for three field orienta- 
tions and for three angles 9 between field and current direc- 
tion. The parabola is a fit of the expression for the diamag- 
netic shift of the B"*" level to the data points. The data points 
are derived from curves as in the inset, which shows the field 
dependence of the resonance of device 1/888/1 (the field step 
is 2 T). 



remove this hole from the ion to the valence band edge, 
is about 2.0 meV. This is close to 0.055Ry* = 2.5 meV 
predicted by the ion model, which is successfully used 
for the D~ state^. Here Ry* is the effective Rydberg for 
boron in silicon, which equals the ground state energy 
El = 45.7 meV^. Since the measured resonance comes 
from the (5— layer, equal parts of the bias drop across the 
barriers at either side of the Coulombic potential well 
associated with a boron impurity. Thus, the resonance 
voltage is K-es = 2[<j)B- E+{Q)]/e [see Fig.lHb)]. We de- 
termine (pB from the temperature dependence of the zero- 
bias resistance Rq. When plotted versus l/T, the loga- 
rithm In (T^Rq) clearly shows activated behaviour in the 
range 15-20 K. Interpreting this as Richardson-Dushman 
thermionic emission^ ■'■ of holes over the barrier, we find 
(j)B — 11.7 meV from a fit to the data. This is not too far 
from (f>B — Ai?i, — Ep = 8 meV, which we obtain from 
the estimatoi^ for AEy{NB) and from photoluminescence 
measurementai^ of Ep{Nb), each subject to uncertainty. 
The values (/-b = 11.7 meV and E+{{)) = 2.0 meV predict 
Vrcs — 19.4 mV. This deviates strongly from the mea- 
sured position and thus seems to exclude the B+ state. 
However, the concentration Nb = 5 x 10^"^ cm~'^ of the 
(5— layer is high enough for the tail of the wave function 
of the second hole to be appreciable at the nearest B" 
atomsi^, so that the B+ ions are not isolated. The addi- 
tional Coulomb attraction of these B" atoms and spread- 
ing of the electron charge among them (reducing the hole- 
hole repulsion at the B+ state) cause a stronger bind- 
ing. This effect increases with increasing concentration. 
The measured Kcs ~ 10 mV gives E'^{Q) w 6.7 meV. 
As demonstrated in Fig. 01 this value nicely falls in the 
bandwidth obtained by extrapolating the scarce exper- 
imental datai^4i^*i& on the concentration dependence of 
the ionization energy for B"*" and D~. 
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Figure 4: Plot of experimental ionization energies for B"*" 
and D~ versus doping concentration (from literature). Above 
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10 cm the ionization energy increases. The data 



point of the present work (♦) nicely falls in the extrapolated 
bandwidth which indicates the trend of the literature data. 



In the lower part of the temperature range of Fig. 2 
the resonance is much wider than the theoretical width 
3.5fcT of a Fermi-smcared sharp resonance, implying a 
rather large zero temperature width (w 1.5 meV). Mecha- 
nisms contributing to this are life-time broadening, disor- 
der broadening and broadening due to the finite width of 
the (5— layer. Life-time broadening determines the intrin- 
sic width of a resonance coming from a single impurity. 
For the B^ level life-time broadening is appreciable, since 
it is so close to the ionization level. The effect is enhanced 
by the large bias field at resonance (2.5 kV/cm), which 
weakens the collector barrier and thus shortens the life 
time. We will come back to this in discussing broadening 
of the resonance with increasing magnetic field (see be- 
low). Disorder broadening arises from the different local 
surrounding of B*' centers by adjacent B*^ centers, which 
causes a distribution of levels. The finite width of the 
(5— layer implies a range of values of Vres, since the posi- 
tion of the atom defines the barrier thicknesses and thus 
the distribution of the bias over them. Without correct- 
ing for a decay of the resonance for off-center positionsii, 
a range for V^cs of 3.5 mV is derived from only the finite 
(5— layer width. We take this as a sign that the i5— layer 
width is an important source of broadening. 

The resonance shift to higher bias (Fig. reflects the 
diamagnetic shift Ai?dia(S) of the B"*" level in a mag- 
netic field. This shift, termed diamagnetic because of 
the related negative susceptibility, is towards the valence 
band edge, in agreement with the observed peak shift 
to higher biases. In magneto-optical spectroscopy the 
diamagnetic shift is not obtained directly, since it is ab- 
sorbed in the field-dependent binding energy E~^(B) = 
E+{Q) + ^huj ~ A£'dia(B). Here ifiw is the energy of 
the first Landau level, which is the valence-band edge in 
field. In our experiment, however, since the emitter dop- 
ing level Nb = 10^^ cm~^ is high enough to block Landau 
level formation in this layer, the emitter Fermi level is a 
field-independent reference energy enabling direct mea- 
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surement of AEdia.{B). For the weak fields used here, 
first order perturbation theory estimates the diamagnetic 
shift asi^ AEdUB) = ^ Y^^jf- Here m is the effec- 
tive mass of the holes of the B+ ion and is the mean 
square distance of the z-th hole to the B~ core {i = 1,2). 
TakingiS mih — O.lSmo and to/j/i — 0.54mo for the light 
and heavy hole hole masses, respectively, and r\ ~ Cg for 
the first hole (ao=3.9 nm is the Bohr radius of the 
atom), we fit the expression for A£'dia(i?) to the complete 
set of data points of Fig. O This yields the fit shown in 
the figure and corresponding values r| = (1.6ao)^ and 
r| = (3.4ao)^. The range defined by these values is con- 
sistent with r2 = 2.4ao cited for the B+ stateSS. We 
attribute the deviation of the fitted curve from the ex- 
perimental trend to the use of an atomic physics model 
for Ai^dia in a sohd state system. 

Broadening of the resonance with increasing magnetic 
field may be unresolved splitting and/or life-time broad- 
ening induced by the Stark effect. Since the resonance 
voltage increases with increasing magnetic field, the elec- 
tic field at resonance, and thus the Stark broadening, 
increase as well. Data for Stark broadening of the _B+ 
level are not available. Therefore, we take as a measure 
the broadening of the far infrared absorption line due 
to transition from the ground state of B^ atom in sil- 
icon to the first excited state, which was measured up 
to 1.0 kV/cm2i. Extrapolation of the data of Ref. [21] 
to the fields of our experiments yields an increase of 
the halfwidth between T and 14 T of 0.1 meV, to be 
compared with our measured increase of 0.5 meV. Stark 
broadening apparently plays a role. 

Finally, we discuss the background contribution to the 
conductance (see Fig. 2), which shows a weak tempera- 
ture dependence below 4.2 K. The parabolic shape of the 
background suggests direct tunneling as transport mech- 
anism. However, the conductance at = is several 
orders of magnitude higher than expected for direct tun- 



nelingS^, so that it is excluded. Hence, the background 
conductance is due to hopping resulting from the back- 
ground doping in the barrier. This hopping is thermally 
activated close to = and field activated at higher 
biases {\V\ > 2 mV). For biases exceeding the barrier 
height {\V\ > 11.7 mV), the barrier becomes increasingly 
weak for hopping and finally for direct tunneling, giving 
a further conductance increase. Above 4.2 K the conduc- 
tance in the range {\V\ < 2 mV) undergoes a transition 
to coexistence of thermal hopping and thermal activation 
over the barrier and finally to dominance of activation 
over the barrier. 

In conclusion, we have studied resonant tunneling 
through a Si barrier 5— doped with boron impurities. 
The conductance resonance observed is due to tunneling 
through the B+ state of the impurities. The structure of 
our device enables direct observation of the diamagnetic 
shift of the B+ state. The measured magnitude of the 
shift agrees well with the theoretical description, yielding 
a proper value of the orbit size for the second hole of the 
B+ state. The binding energy of the B+ state turns out 
to be enhanced as a result of overlap of the wave function 
of the second hole of the B'^ state with the nearest boron 
impurities. Our next step will be miniaturizing the de- 
vices to the level of one dopant atom (diamcter«50 nm), 
enabling studies of the effect of wave-function manipula- 
tion on transport through the atom. 
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